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1.  INTRODUCTION 

The  investigation  of  both  theoretical  effects  and  computer  simulation  models  of  the  collection  of 
ion  current  by  a  Langmuir  probe  in  the  plasma  wake  of  a  shadowing  body  is  the  focus  of  this 
report  Results  of  such  experimentation  have  been  used  in  flight  predictions  [Cooke  et.  al. ,  1994] 
of  the  WSF-CHAWS  STS-60  spaceflight  mission.  It  has  been  found  that  the  collection  of  ion- 
current  in  a  plasma  wake  is  a  space-charge-limited  process  with  a  heavy  dependence  upon  an 
angular-momentum  mechanism  Previous  attempts  [Enloee/.  al 1993,  Cooke  et.  al 1994]  have 
been  made  to  model  effects  of  scaling  plasma  parameters  and  pose  theoretical  models  to  asses  the 
variations  of  current  collection  with  the  variations  of  plasma  parameters  Although  expected  by 
orbit-limited  models  of  ion  collection  in  the  plasma  wake,  the  magnitude  of  the  collected  current 
is  not  a  linear  function  of  plasma  temperature  (T),  plasma  density  (N),  nor  ion  Mach  number  (MV 

2.  THEORETICAL  MODELS 

Given  that  the  angular  momentum  of  the  particles  is  a  dominant  characteristic  of  particle 
collection  in  the  wake  [Enloe  et.  al ,  1993],  we  must  speculate  on  the  impact  of  the  scaling  of 
particular  plasma  parameters  for  both  orbit-limited  and  space-chargc-limited  analytical 
approaches  For  the  scope  of  our  study,  we  look  at  the  first  order  impact  of  the  most  dominant 
plasma  parameters  plasma  density,  temperature,  Mach  number,  and  atomic  mass  number  In  any 
plasma  environment,  we  know  that  the  ion  thermal  velocity  is  proportional  to  the  square  root  of 
the  plasma  temperature  to  mass  (AMU)  ratio,  more  conveniently  written 


vth 


(kT/mi)l/2  x  (T/AMU)1/2 


(2  1) 


For  a  moving  object  in  a  stationary  plasma  or  a  stationary  object  in  streaming  plasma,  we  know 
that  the  relative  plasma  stream  velocity  (vq),  ram  current  density  Gram)-  antl  plasma  beam 
energy  (Eq)  are  related  as 
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where  M  is  the  mach  numbei  of  the  streaming  plasma  relative  to  the  object  immersed  in  the 
plasma  From  orbit-iimited  theory,  we  would  expect  that  an  increase  in  the  ram  current  density  by 
a  factor  of  two  would  be  followed  by  an  increase  in  current  collection  by  a  facor  of  two  From 
space-charge-limited  theory  for  a  body,  such  as  the  WSF-CHAWS  experiment,  where  the 
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dimensionless  characteristic  radius  (r=RA[))  is  much  greater  than  unity,  the  current  collection  of  a 
probe  in  a  plasma  is  limited  by  the  Debye  shielding  of  the  plasma  sheath  set  up  by  the  probe 
potential  [Langmuir  and  Blodgett,  1924],  Otherwise,  we  can  assume  that  the  current  collection  in 
the  wake  of  the  WSF-CHAWS  experiment  is  orbited  limited  when  the  dimensionless  characteristic 
radius  is  less  than  unity.  Such  a  case  can  occur  when  the  density  is  at  a  critical  minimum,  thus 
making  the  influence  of  the  Debye  length  more  appreciable  to  the  characteristic  radius  (R)  of  the 
experiment  In  the  case  of  the  WSF-CHAWS  experiment,  a  rani  density  on  the  order  of  104 
becomes  highly  unlikely  for  a  body  in  LEO  Therefore,  we  can  assume  that  a  simple  scaling  of 
the  ram  flux  density  and  the  associated  collected  wake  current, as  in  accordance  with  the  previous 
orbit-limited  assumption,  once  again  is  inaccurate 


2.1.  Mach  Number  and  Plasma  Temperature  Variations 

Reverting  to  a  more  basic  approach,  we  must  first  concentrate  on  known  relations  In  congruence 
with  the  WSF-CHAWS  spaceflight  operation,  we  know  that  the  orbital  speed  of  the  shuttle  and 
WSF  is  nearly  constant  This  constant  velocity,  for  any  given  species  of  plasma,  denotes  that  the 
streaming  velocity  (beam  energy)  of  the  plasma  also  will  be  a  constant,  relating  T  and  M 
appropriately  by  Eqn  2  4 


M2*-  (2  5) 

To  model  the  phvsical  effects  of  the  variations  of  these  parameters,  we  look  at  the  ram  flux  from 
the  origin  of  the  plasma  wake  disturbance  at  the  edge  of  the  WSF-CHAWS  shadowing  body 
Since  the  WSF-CHAWS  experiment  is  traveling  in  the  -Z  direction  at  a  constant  velocity  vzo  and 
a  transverse  ion  expansion  into  the  wake  with  a  mean  velocity  of  vxo,  we  have  the  corresponding 
distribution  function 
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Integrating  this  function  with  respect  to  each  velocity  component 
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which  agrees  with  Parker  [  1980]  Analysis  of  this  result  provides  important  information  about  the 
effect  of  the  plasma  temperature  (Mach  number)  on  the  ram  current  density  For  Mach  numbers 


2 


above  the  value  of  two  in  Eqn  2.8,  the  error  function  saturates  to  unity  and  the  negative 
exponential  saturates  to  null,  indicating  that  the  flux  distribution  is  unaffected  by  Mach  numbers 
above  the  value  of  two.  In  comparison  to  oxygen  ions  in  LEO,  the  WSF-CHAWS  experiment 
operates  in  a  Mach  number  range  above  the  value  of  seven.  Therefore,  the  collection  of  current 
due  to  oxygen  ions  on  a  probe  in  the  WSF-CHAWS  wake  is  invariant  to  perturbations  in  the 
plasma  temperature  For  hydrogen  ions  in  the  same  environment,  the  WSF-CHAWS  experiment 
operates  in  a  Mach  number  range  near  the  value  of  two,  indicating  that  the  ram  density 
characteristics  may  be  dependent  upon  the  plasma  tempeiature 

2.2.  Plasma  Density  Variations 

Although  variations  in  the  plasma  temperature  (Mach  number)  have  little  effect  on  the  I-V 
char  icteristics  of  the  Langmuir  probe  in  the  WSF-CH  AWS  wake,  variations  of  the  ambient  ion 
density  of  a  LEO  environment  exist  As  Eqn  2  8  shows,  the  flux  density  is  directly  proportional 
to  the  density  of  the  plasma  As  stated  before  by  the  orbit-limited  theory,  we  can  assume  that  the 
current  collection  in  the  wake  will  be  proportional  to  this  density,  but  first  we  must  look  at  the 
effect  of  the  density  on  the  space-charge-limited  collection  of  a  Langmuir  probe  according  to 
Langmuir  and  Blodgett  [1924]  From  the  well-known  Child-Langmuir  law  of  space-charge- 
limited  current  between  charged  concentric  spheres,  we  have 
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where  eG  is  the  permittivity  of  free  space,  V  is  the  potential  difference  between  the  spheres,  and  D 
is  function  of  the  radii  of  the  spheres  Using  the  three  dimensional  flux  density  to  a  surface  (let  J0 
-  eN[kT/27tM]  Debye  length,  and  dimensionless  potential.  Eqn  2  9  can  be  written 
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From  the  Park  ir  [  1980]  tabulation  of  Table  II  of  Langmuir  and  Blodgett  ( 1 924], 
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where  H0  is  the  characteristic  radius  of  the  Langmuir  probe  and  Rs  is  the  sheath  radius 
Incorporating  Eqn  2  10  into  Eqn  2  I  I ,  yields 


(2.12) 
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Using  2. 12  as  the  radius  of  the  spherical  sheath, 
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From  the  definition  of  the  ion  Debye  length  and  the  proportion  of  Eqn.  2.1,  the  current  to  the 
Langmuir  probe  is  proportionally  related 


I  =  4nJ„Rs2 
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Therefore,  we  see  that  for  a  non-flowing  plasma,  the  space-charge-limited  current  varies  with  the 
3/7  ths  power  of  the  variation  of  the  plasma  density  The  -^/ j  4  ths  power  variation  associated 
with  the  plasma  temperature  is  considered  ineffective  via  the  derivation  of  Eqn  2  8  To  apply  a 
similar  reasoning  to  the  plasma  stream  that  evolves  about  the  edges  of  the  shadowing  body  proves 
to  be  rather  difficult  We  cannot  assume  a  spherical  sheath  as  in  Eqn  2  13,  nor  is  the  plasma 
density  in  the  wake  uniformly  distributed  about  a  spherical  probe  What  is  important  is  the  power 
relation  of  the  density  to  the  Langmuir  probe  current  On  the  other  hand,  if  we  assume  that  the 
collection  potential  of  a  Langmuir  probe  in  a  plasma  wake  is  sufficiently  high  enough  (V  »  kT/e) 
above  some  collection  threshold,  we  can  (for  the  moment)  disregard  the  angular  momentum 
properties  of  ion  trajectories  in  the  sheath  and  use  the  Langmuir-Blodgett  relation  tor  the  current 
scaling  with  density  Such  a  situation  is  physically  unrealizable  because  the  sheath  edge  is  not 
well  defined  and  the  initial  velocities  of  the  ion  trajectories  into  the  sheath  are  highly  dependent 
upon  the  angular  momentum  of  the  ion  panicle  From  this,  we  can  conclude  that  the  current 
collection  mechanism  in  the  plasma  wake  is  neither  orbit-limited  or  space-charge-limited  Rather, 
it  is  in  a  quasi-space-charge-orbit-limited  regime  where  the  current  scaling  via  the  density  can  be 
described  mostly  by  the  relation  of  Eqn  2  14  in  the  space  charge  sense,  but  the  particles  also 
exhibit  some  behavior  relevant  to  angular  momentum  concepts 
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3.  EXPERIMENTAL  RESULTS 

We  have  used  the  POLAR  [Lilley  et.  al.,  1985]  computer  simulation  to  model  the  chamber 
experimentation  [Enloe  et.  al.,  1993]  and  the  WSF-CHAWS  STS-60  mission  [Cooke  et.  al., 
1993]  with  good  results.  The  POLAR  code  was  able  to  model  the  chamber  experiment  as  shown 
in  Figure  1.  The  in  situ  results  [Enloe  et.  al.,  1993]  and  POLAR  simulations  are  shown 
respectively  in  Figures  2  and  3  The  identification  of  the  "sweet  spot"  of  orbiting  ions  and 
similarities  in  the  magnitude  of  current  collection  provided  a  reasonable  validation  A  similar 
geometry,  shown  in  Figure  4,  was  used  to  model  the  parametric  variations  more  closely  by 
removing  the  effects  of  the  ion  circulation  about  a  Langmuir  probe  with  quasi-radial  symmetry. 
The  small  effects  of  Mach  number  and  plasma  temperature  variations  are  shown  in  Figure  5.  The 
effects  of  plasma  density  variations  and  the  Langmuir-Blodgett  scaling  assumption  is  shown  in 
Figure  6  Similar  relations  of  the  plasma  density  variations  are  also  present  in  the  Cooke  et.  al. 
[1994]  data 

4.  CONCLUSIONS 

From  the  theoretical  assumtions  and  predictions  of  the  POLAR  simulation  models,  we  see  that 
varaitions  in  the  Mach  numbers  and  plasma  temperature  have  little  effect  on  the  ion  current 
characteristics  of  collection  when  the  Mach  number  is  greater  than  two  It  was  found  that  the 
collection  of  ions  to  a  Langmuir  probe  in  the  plasma  wake  of  a  shadowing  body  is  a  space  charge 
limited  process  that  is  heavily  dependent  on  the  angular  momentum  of  the  particles.  The  collected 
ion  current  is  expected  to  scale  with  variations  in  the  plasma  density  to  the  ^/j  ths  power. 
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Figure  1  POLAR  Chamber  Simulation  Shadowing  Body  and  Langmuir  Probe  Configurations: 
a  )  side-view  (x-z)  of  POLAR  object,  b.)  three  dimensional  view 


Figure  2  Current  collected  as  a  function  of  negative  potential  on  the  Langmuir  probe  for  varying 
axial  positions  The  change  in  the  shape  of  the  1-V  curve  for  z  =  3  cm  is  due  to  the  effect  of  the 
ions  orbiting  the  sphere  [Enloe  et.  al 1993] 
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Figure  3.  Current  collected  as  a  function  of  negative  potential  on  the  Langmuir  probe  for  varying 
axial  positions  from  the  POLAR  simulation. 


a)  b) 

Figure  4  POLAR  simulation  shadowing  body  and  Langmuir  probe  configurations  used  to  limit 
effects  of  ion  circulation  a  )  side-view  (x-z)  of  POLAR  object,  b  )  three  dimensional  view 


POLAR  A  (pos.  3)  N-3e+13 


Figure  5.  I-V  collection  characteristics  as  a  function  of  plasma  temperature  and  related  Mach 
number. 


a.)  No  density  scaling  b.)  Langmuir -Blodgett  density  scaling 


Figure  6  I-V  collection  characteristics  as  a  function  of  density.  The  single  element  Langmuir 
probe  was  replaced  by  an  elongated  probe  extending  from  the  shadowing  body  to  prevent 
trajectory  orbits 


